Optical spectral weights and the ferromagnetic transition temperature of CMR 
manganites: relevance of double-exchange to real materials 
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We present a thorough and quantitative comparison of double-exchange models to experimental 
data on the colossal magnetoresistance manganese perovskites. Our results settle a controversy 
by showing that physics beyond double-exchange is important even in LacrSrcaMnOs, which has 
been regarded as a conventional double-exchange system. We show that the crucial quantity for 
comparisons of different calculations to each other and to data is the conduction band kinetic 
energy K, which is insensitive to the details of the band structure and can be experimentally 
determined from optical conductivity measurements. The seemingly complicated dependence of Tc 
on the Ifund's coupling J and carrier concentration n is shown to reflect the variation of K with 
J, n and temperature. We present results for the optical conductivity which allow interpretation 
of experiments and show that a feature previously interpreted in terms of the Hund's coupling was 
misidentified. We also correct minor errors in the phase diagram presented in previous work. 



o 
o 



> 

in ■ 

o : 

m 
oo 

o ■ 

o^ ; 



o 
o 



X 



I. INTRODUCTION 

The colossal magnetoresistance (CMR) rare earth 
manganese perovskites first attracted attention in the 
1950s because of the range of magnetic and structural 
transitions they display. Recent interest has been re- 
vived by the extremely large ("colossal") magnetoresis- 
tance displayed by some, CMR materials, coupled with 
their rich phase diagramEJ. Despite this long and intense 
study, much of the physics of CMR remains controversial, 
with basic issues still subject to debate. In this paper we 
address two such issues: the first is in what sense the 
standard 'double-exchange only' model (defined below) 
describes the physics in the regions of the phase dia- 
gram where the ground state is a ferromagnetic metal. 
The second is the interpretation of the optical conduc- 
tivity spectrum, and in particular which portions of the 
observed spectrum pertain to the conduction band elec- 
trons responsible for the interesting physics of the CMR 
materials. In addition to its direct relevance to CMR, 
we believe our work is of broader significance for the the- 
ory of interacting electrons in solids, as a contribution 
to the fundamental issues of the quantitative comparison 
of many-body calculations to the properties of materi- 
als, and to the interpretation of the optical spectra of 
correlated electron materials. 

It is generally agreed that a crucial aspect of CMR 
physics is 'double-exchange '□: the mobile carriers (which, 
for CMR materials, are Mn Cg symmetry d-electrons) are 
strongly coupled ferromagnetically to localized core spins 
{Mn t2g symmetry electrons). This means that core spin 
alignment dictates carrier motion, which in general leads 
to ferromagnetism. What, if any, additional physics be- 



yond double-ifxchange is important for the manganites is 
still debatedH. For example, some authors have argued 
that double-exchange only models correctly predict the 
magnetic transition temperatures of CMR nraterialsBu, 
while others have argued, that they do notd'u. Pub- 
lished calculations of TcEIQ have not resolved the issue 
because the results depend on model parameters such as 
bandwidth, interaction strength and carrier density, of- 
ten chosen arbitrarily or in a manner inconsistent with 
the physics of the manganites. We attempt to settle 
the question by establishing the full phase diagram and 
by presenting precise calculations of within a specific 
model. Most importantly, by showing how the parame- 
ters on which the calculated Tc depends can be related to 
measured properties of real materials and to parameters 
used in other calculations. 

Another important set of issues concerns the optical 
conductivity, which in the CMR materials has a strong 
dependence on chemical ^anmosition, frequency, temper- 
ature and magnetic fieldEJij E3. It is widely believed that 
information extracted from these data will be helpful in 
elucidating the physics of CMR, but this goal has not yet 
been fully achieved. Even the basic question of which 
parts of the observed spectrum are due to the,cjOiiduc- 
tion band degrees of freedom is still not settledl3~E3. We 
show here that a comparison of the data to the theoreti- 
cally calculated magnitude and temperature dependence 
of optical spectral weight, and of its relationship to Tc 
can resolve this i^Wie. Our results indicate that several 
previous analysesou underestimate the actual conduction 
band spectral weight, and suggest that-a previous paper 
by several authors including one of usE2l misidentified a 
key feature in the data. This misidentification arose pri- 
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marily from reliance on a theoretical calculation based 
on inappropriate system parameters. 

In this paper we perform a systematic and quantitative 
comparison of the predictions of double-exchange-only 
models to data. We incorporate "real materials" aspects 
via a tight-binding fit to the band structure; in these 
materials a simple nearest neighbour hopping model rep- 
resents the real material parameters adequately. In par- 
ticular, we show that the crucial quantity for compari- 
son of calculations to data or to other calculations is the 
conduction-band kinetic energy, a local expectation value 
related to the electron hopping amplitude and defined 
more precisely in Eq. (||). Because this kinetic energy is 
a local quantity, it can be easily and reliably calculated 
and is insensitive to the fine details of the band structure. 
When expressed in terms of the kinetic energy, the appar- 
ently disparate results of a variety of computations are 
seen to be in quantitative agreement. More importantly, 
because the kinetic energy may ha determined from op- 
tical conductivity experimentsOEj, our results enable a 
detailed comparison of theory to experiment. 

To treat the many-body physics of double-exchange, 
we usa— .the highly successful dynamical mean-field 
methodtS. This provides a detailed and apparently reli- 
able treatment of the local correlations which are crucial 
to the conclusions we wish to draw. We will argue that 
the effect of the omitted physics, which has to do with 
intersite fluctuations, may be reliably estimated for the 
issues of interest to us and is not too large. 

The application of the dynamical mean-field method 
to double-exchange systems was pioneered by FurukawaO, 
whose important work established the basic usefulness of 
this method and presented many valuable results con- 
cerning the phase diagram, conductivity and other prop- 
erties. However, his work was incomplete in some re- 
spects, incorrect in some details, and in many cases em- 
ployed physically irrelevant parameters. The full treat- 
ment we present here is therefore needed. 

We believe that the analysis presented in this paper 
may be useful in the general context of the development 
of a predictive theory of dynamical and ordering phe- 
nomena in systems, such as transition metal oxides, with 
strong interactions. This is an important goal of elec- 
tronic condensed matter physics; and recent advances in 
computational power and in the techniques of many-body 
physics suggest that it-^nay be obtainable. The dynami- 
cal mean field methodt^ seems particularly promising in 
this respect. It is computationally tractable, includes in- 
coherent and inelastic effects, and can be combined with 
conventional band theory. Investigations of the extent to 
which this method is useful in the computation of real 
materials properties are urgently needed. The work de- 
scribed below is such an investigation. In this context we 
point out that one of the great advances in the theory 
of equilibrium critical phenomena arose from a careful 
and detailed application of mean- field theories. We be- 
lieve that a similar opportunity may be present in the 
application of the dynamical mean field theory to strong 



correlation physics. 

The rest of the paper is organized as follows. Section 
II defines the double-exchange-only model, and explains 
how its parameters should be defined and related to ex- 
periments. Section III presents the dynamical mean-field 
formalism. Section IV presents the numerical and ana- 
lytical results for the phase diagram, and corrects what 
seems to be a minor error in the phase diagram proposed 
inlll. Section V presents results for the spectral weight 
(integrated area) in different regions of the optical con- 
ductivity spectrum. Section VI provides a summary of 
our theoretical results and a detailed interpretation of ex- 
perimental data, particularly of optical conductivity and 
its relation to Tc- Readers interested only in this aspect 
may proceed directly to this section. Section VII is a con- 
clusion and presents some possible directions for future 
research. 



II. MODEL 

The double-exchange-only (DE) model involves elec- 
trons (with orbital indices a, b) moving in a band struc- 
ture defined by a hopping matrix tfj and a chemical po- 
tential /Lt, and connected by a Hund's coupling (J > 0) 
to core spins S, which we take to be classical. We denote 
the operator creating an electron of spin a on orbital 
a of site i by dl^^ arid define the double-exchange only 
Hamiltonian H]je byll3 

HdE = - X! ^i]4aad3ba ~ fJ'^dl^^diaa 

<ij>,ab iaoL 

-J^S.,- dl^^aapd,a/3 ■ (1) 

iaaf3 

Here t^j is the amplitude to hop from site i, orbita a 
to site j, orbital b. The calculated band structureEll is 
well fit by a tfj which involves only nearest-neighbour 
hopping. 

A crucial quantity is the electron kinetic energy K, 
defined by 
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where Cp** is the dispersion obtained by Fourier trans- 
forming tfj, Z is the number of nearest neighbours and 
d is the spatial dimension. In the limit J — > cxd, K is 
the only relevant energy scale in the model and in par- 
ticular the magnetic transiti|Qu-.t|eHiperature Tc has been 
found to depend only on ida'tsO. K is also a suitable 
quantity for comparison to experiment because, in the 
physically realistic limit of nearest-neighbour hopping, 
K may be determined from an analysis of the optical 
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conductivityll3jlj. We would like to mention that this K 
includes a sum over all bond directions and is 3 times 
larger than the K defined in Quijada et al. 

III. FORMALISM 

In this paper we study the magnetic transition tem- 
perature and kinetic energy of Hjj e as functions of car- 
rier density n and Hund's coupling using the dynamical 
mean field approximation, supplemented by an analyti- 
cal treatment of the Jj-,-^ limit and by cornparison to 
results of Monte CarldZl and series expansionsEj. 

The dynamical mepi field method has been discussed 
extensively elsewhereE-3 and therefore we omit formal de- 
tails. This scheme was first obtained in the limit d —^ oo 
with tij ^ l/Vd, and is often referred to as the infinite 
dimensional limit, but the central approximation is ac- 
tually the neglect of the momentum dependence of the 
electron self-energy, E. While this becomes exact in tlie 
limit d ^ oo, it is quite a good approximation in d = 3li3, 
where the momentum dependence of E is quite weak. 

To make the importance of the neglect of momentum 
dependence clear and to establish notation we sketch 
a derivation of the dynamical mean field equations. A 
momentum independent S implies that all of the many- 
body physics may be derived from the local (momentum- 
integrated) Green's function Gioc, defined by 



(3) 



where the Green's function, a matrix in orbital [ah) and 
spin {aj3) space, is 



G{p,u:)=[u: + ^,-ef~Y.tp{u^)\ ' 



(4) 



It is convenient to define the density of states summed 
over bands N{e) — J2\ J '^p)^ where are the 

eigenvalues of e"''. We shall be interested in the pseudo- 
cubic manganese perovskites in the doping regimes where 
orbital order is not important. For these materials 
Nx{e) = / j§^S{e ~ Ep) is independent of A; thus we 
use N(e) = norbN\{e), where we have introduced the 
number of orbitals Uorb- For the physically relevant case, 
Uorb = 2, but-jthe case Uorb = 1 has been considered by 
other workeraj, so we retain Uorb as a variable to facilitate 
comparisons. The local Green's function is then propor- 
tional to the unit matrix in ab space, G = Gioc.affSab and 
the coefficient is 



deN{e)- 



1 



Sq/3(w) 



(5) 



Gioc depends only on frequency and is therefore the so- 
lution of a single-site problem. This problem is specified 
by a mean-field function b^(uj) whichj-4p the model of 
present interest, is related to Gioc viaS'tZl 



Gloc,ap{^) — / 

with 



p{n) 



-,n^s/ s (6) 



Pin) = exp 



bal3i^) — JS ■ (Ja0 

Tr In bap (iujn ) - JS ■ a^p 



(7) 



and Zioc ensures that / d^VlP{n) = 1. The self-energy 
is defined by 



T^apiu) = bap{uj) - Gi^\^^p{uj) 



(8) 



and the mean- field function bap{ijj) is fixed by the re- 
quirement that Eq. (H) holds with Gap{p,(^) defined by 
Eq. d), E defined by Eq. §) and Gioc by Eq. §). 
Using the momentum independence of the self-energies 
within this mean-field theory, the electron kinetic energy 
(Eq. (0)) can be written as 



z ^ 



^—^epTr[G{j>,iu„ 



jTj^Tr [G 

IOC 



(9) 



In a ferromagnetic state with magnetization direction 
m we have 



hap{i~o) = bo{uj) + 6i(w)to • aap 



(10) 



{bi{uj) = in the paramagnetic state). If the spin axis is 
chosen parallel to m then bap becomes a diagonal matrix 
with components parallel (&| = 5o + 6i) and antiparallel 
(fii = bo — bi) to m. 

The precise form of the equation for b^ depends on 
the form of the density of states A''(e), but the important 
behaviour of the physical observables do not, as long as 
N{e) has a finite second moment (/ deN{e)e'^ < oo). We 
perform our calculations on the Bethe lattice with a semi- 
circular density of states per orbital 



N{e) = 



27r<2 



(11) 



Here t is a bandwidth parameter; the full bandwidth of 
the non- interacting (J = 0) problem is 4i. This cor- 
responds to number of nearest neighbours Z — 2d, as 
Z oo with tZ fixed. For this N{e), the self-consistent 
equations for the feo- are 



6| = iuj„ + + J d{cose)P{0) X 

{-bi + JcosO) 
b^bi- {b^ -bi)Jcose- J"^ 



bi — iuJn + - 



d{cos9)P{e) X 



(6| + JcosO) 



b]bi - (6| - bi) JcosO - J2 



(12) 
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and the local Green's function (Eq. ||) can be written in 
terms of the mean-fields as 



(13) 



The self energies are evaluated using Eq. (g) and the full 
Green's function from Eq. (^. 

We will also be interested in the optical conductivity, 
a. The contributions to a from the states described by 
H]j E are obtained by couplng an electromagnetic field to 
Hde (Eq- 0) via the Peierls coupling tij ij^e*'^^''*"''^^ 
and performing linear response. If the self-energy is mo- 
menjttim independent then there are no vertex corrections 
in ac3 and 



(14) 



where jp is the current operator determined by the 
Peierls substitution; it includes both intra and inter- 
band transitions within the manifold of conduction band 
states. Since we concentrate on materials with cubic sym- 
metry, (7ij ^ Sij and hereafter we suppress the spatial 
indices on a. 

This conductivity obeys the sum ruleP"^ 



S'(oo) 



e2 







n 



iSaba 



H.c. 



(15) 



If the nearest neighbour hopping is dominant then 
/q°° § dri(j{^l) = -^^rrjK. The experimentally measured 
optical conductivity includes additional transitions not 
described by Hue and obeys the familiar f-sum rule 
^ dVla{Vl) = with n the total (conduction and 
valence) electron density and m the bare mass. Extract- 
ing the portion pertaining to the Cg electrons from the 
measured conductivity has been controversial. We will 
show below how to do this using the information on spec- 
tral weights we derive below. Here we concentrate on 
spectral weight defined in Eq. ( p^ ) . The real part of the 
optical conductivity (Eq. |lj) can be written as 



a{n) 



/ defciV(efe)^(efe) 



duj f{u)- J{Vl + u}) 



X A^(efc,w)A^(efe,r2 + w) 



(16) 



where A^{ek,uj) = —-^ImGaisk,^) is the spectral func- 
tion, N{ek) the density of states (Eq. ^ and (p{ek) — 
(4^2 — e|)/3 is the current vertex for the Bethe lattice. 
We obtain (p{ek), the form of whicii,for the Bethe lat- 
tice has not been given before (Ref.tj presents (p for the 
Gaussian density of states) , by requiring that K be given 



both by the explicit integration of a (Eq. 16) and by the 
generally valid expression 

/ dek—f{uj)ekN{ek)ImG{ek,uj). (17) 



Requiring the two to be equal yields the differential equa- 
tion 



d 

- — [N{ekMek)] - N{ek)ek , 
the solution of which yields '^{tk)- 



(18) 



IV. NUMERICAL METHODS 

We now outline the numerical methods used to solve 
Eqs. (|2|) and compute physical quantities. For com- 
putational convenience, we rewrite the self-consistency 
Eqs. ( p^ ) in terms of 6o and bi defined in Eq. ^ At 
T > Tc, there is no magnetic order and bi{uj) — 0. At 
T — and non-zero J, all spins are ferromagnetically 
aligned implying P{9) — > d{cos9—l), and one can analyt- 
ically solve for the bo,i{i(^n) from Eq. (|l2|). We initialize 
feo.i with the T = expressions, and solve Eq. ( p2|) on 
the Matsubara axis by direct iteration. 
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FIG. 1. The dependence of Tc on the Matsubara cutoff 



for n — 0.75, Uorb ~ 1 and J — lOt; 



ranges from 



(4t -|- J) to 4(4f + J). The quadratic dependence on t/ujmax 
holds for the well. For our calculations we choose 

^max = 3(4f -f J) which produces a value of Tc that is within 
0.5% of the LUmax oo value. 

It is convenient to compute P{0) by solving Eq. ( [l^ ) on 
the Matsubara points, and then use this converged P{0) 
to solve the equations on the real axis. Convergence to 
within an rms error of 10~^ was usually achieved within 
12 iterations unless one is close to the magnetic Tc where 
there is critical slowing down. The electron density n and 
the normalized magnetization density m are given by 
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m= I d{cose)coseP{6) . 



(19) 
(20) 



An accurate value of the transition temperature Tc was 
most conveniently obtained by computing values of m 
(Eq. 1^) in the range 0.17 < m < 0.3 and finding by 
fitting to the mean- field expression im?{T) = a{Tc — T), 
with a and fit parameters. 

The next issue is that of the number of Matsubara 
points needed to compute Sioc{ff), n and K. At large ujn 
, the asymptotics of the mean fields are 



bo = iujn + M - 



1 + 



iuJn + ^ {iUn + m)^ 



hi 



Jm 



{iuJn + M)^ 



(21) 



In our computations we choose a frequency cutoff LUmax 
and evaluate higher frequency contributions analytically 
using the asymptotic form given in Eq. (|^). The errors 
in physical quantities are of order w"^^.. This is illus- 
trated for Tc in Fig. 0. To achieve an accuracy of 5 x lO"** 
in n, it was required to choose ujmax > 3(4i -I- J). Choice 
of this value of LOmax also ensures that the value of Tc is 
within 0.5% of the ujmax — ^ oo value. 
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A. Analytical results 

We first present mean field results and establish the 
phase diagram, then show how to incorporate fluctua- 
tion corrections. We begin with analytical results for 
small ferromagnetic J. At T = the core spins are mag- 
netically ordered with a characteristic wavevector q. The 
polarized core spins produce an effective magnetic field 
on the conduction electrons which polarizes them, lead- 
ing to a magnetization rriq and to a change in energy, 
which for small J is 

SE^^ ™, • m_, + i J (^Sg ■ m_, + H.c>j (22) 

where Xq is the magnetic susceptibility for nonintcracting 
electrons at wavevector q. Minimization of Eq. (^) leads 
to friq = XqSq and SE = —^J^SqXq- Thus at small J 
the system will order at the wavevector which maximizes 

Xq- 

The wavevector at which x is maximal depends on 
band filling and dimensionality. Here, far|-.consistency 
and to facilitate comparison to other wormu we restrict 
ourselves to infinite-d. For very large d, xil) is inde- 
pendent of q except for regions of width 0{1/Vd) about 
q = and the commensurate antiferromagnetic vector 



g = Q = (tt, TT, ...); therefore at d = oo we need to con- 
sidet^X<;=o 1 Xg=Q and the susceptibility at a typical 9, 
O. Fig. g shows the ferromagnetic (bold line), anti- 



Xi 



ferromagnetic (solid line) and local (dashed line) suscep- 
tibilities for J = calculated for a semicircular density 
of states. The phase corresponding to the maximal x 
is the small J ground state. For < n/uorb < 0.195, 
x(g = 0) > XIocXQj ttiis is the range of dopings where 
the DE model has a low J ferromagnetic ground state. 
At small J and intermediate n (0.195 < n/Uorb < 0.35), 
Xioc is largest, implying order at wavevector different 
from both and Q. At d = 00 all such wavevectors 
are degenerate, implying some sort of highly degenerate 
ground state. Finite d corrections will select a particular 
wavevector. We therefore identify the phase as incom- 
mensurate (IC). For n/riorb > 0.35, XQ is largest and the 
small J phase is a commensurate antiferromagnet. 




FIG. 2. The susceptibilites x(9 = 0) (bold line), xiocai 
(dashed line) and xil ~ (''"1 ■•••)) ^ 8' function of the filling. 
The J = ground state is ferromagnetic where x(? = 0) is 
largest, incommensurate where Xiocai is largest, and commen- 
surate antiferromagnetic where XQ is largest. 

Using mean field theory, we may also obtain an ex- 
pression for Tc- For finite range interactions, mean field 
theory is strictly correct only in the limit d — > cxd. We 
will consider finite d corrections below. Writing Eq. ( [2^ ) 
in real space, focussing on a single site i, integrating 
out the conduction electrons and making a mean field 
ansatz yields a self-consistent equation for the polariza- 
tion on site i. For example, for a ferro- or commensurate 
antiferro-magnetic ordering, the effective field at site i, 
heff,i = JJ2i^j Xi~j {Sj), has magnitude independent of 
i and we find 



= -j^ [XQ - Xloc] 



(23) 



An explicit analytic expression for ferromagnetic Tc can 
be obtained by working out the susceptibilities in the 
mean- field Eq. 03 for Q = 0, and one finds that 



5 



J2 



(24) 



We observe Tc ^ at /i = 1 which corresponds to n = 
Tic — 0.195; for n > ric, the small J ground state becomes 
IC. 

Next we consider the case of very large J, focussing 
first on the ground state energy. For J > J*{n), a 
ferromagnetic state would be fully spin polarized, and 
the ground state energy is simply that of the appro- 
priate density of spin polarized electrons moving in the 
relevant band structure, and is therefore of order t. If 
n/uorb = 1 the band is completely full, and the ground 
state energy vanishes, because no hopping is possible. 
For the semicircular density of states, J*{norb) — 2t and 

if n = Uorbi^ - 



the the ground state energy EpM is 



Efm = ~2xt + O {x^t) . 



(25) 



At n = riorb, s. commensurate antiferromagnetic state 
is favoured, because virtual hops between the up and 
down sublattices leading to an energy gain ^ —t^jj are 
possible. For the Bethe lattice with Z nearest neighbours 
and near neighbour hopping i^, the energy is 



Eaf = 



2Zt| 
J 



2<2 



(26) 



where the arrow indicates the o? — > oo limit. Further, 
one may consider an incommensurate state which for the 
purpose of d = cxD energetics has a random spin arrange- 
ment, leading to 



J/t 



t 

Eic = -j- + O {x^t) 



(27) 



1 



PS \ 






FM 


\ic 




Commensurate \ 




AF \ 




0.35 


0.195 



n/n 



orb 



FIG. 3. The phase diagram of the double exchange model, 
as deduced from the analytic arguments at small and large J. 
At small J, one can have ferromagnetic (FM) , commensurate 
anti-ferromagnetic (AF) or incommensurate (IC) phase. At 
large J, the IC phase is not energetically favoured and PS 
indicates phase separation between AF and FM. 



Equating EpM to Eic implies a ferromagnet- 
incommensurate transition at xp-i = ((2 — v^)"/) 
; equating Eaf to Eic yields incommensurate- 
antiferromagnetic transition at xi-a — t/{'sf^J)- Thus 
at small 1 — n the sequence of spatially homogeneous 
phases as was also found at small J and small n. How- 
ever, at small 1 — n and large J, inhomogeneous phases 
are favoured: the standard Maxwell construction applied 
to Eqs. (^,^,^) shows that in this model, which ne- 
glects long ranged Coulomb interactions, the whole small 
X and large J regime is in fact phase separated into com- 
meiisurate AF and F with x ^ 1/2, as previously noted 
byQ. Combining the small J and large J results, we ob- 
tain the phase diagram shown in Fig. (||). The qualitative 
behaviour is correct; we do not have precise numerical 
values for the IC-AF or PS phase boundaries. 

We emphasize that the large- J antiferromagnetism and 
phase separated behaviour is tied to the regime n « riorb, 
where the conduction band of the ferromagnetic state is 
almost completely full. The physically relevant regime 
for the manganites is n < no,.;,/2, where the large J 
behavieur is simply ferromagnetic. Several published 
papersH have asserted that behaviour associated with the 
regime n « Uorb is relevant to the manganites. In our 
view, these assertions are unjustified because Uorb — 2, 
so the physical density corresponds to n/uorb < 0.5. 

It is possible that additional interactions, neglected 
here, could change the effective orbital degeneracy from 
2 to 1. We do not believe this happens in the mangan- 
ites. In support of our view we cite the case of LaMnOa. 
This is an insulator which has a frozen Jahn- Teller dis- 
tortion which acts to quench the orbital degrees of free- 
dom. This material is a (0,0, tt) antiferromagnet: its 
magnetic structure is ferromagnetic planes antiferromag- 
netically coupled. The ferromagnetic in-plane coupling 
means that virtual hopping to the empty orbital is the 
dominant process. The much weaker antiferromagnetic 
bond perpendicular to the plane is believed to be due 
to an additional t2g — i2g superexchange which becomes 
important because the Jahn- Teller order strongly sup- 
presses the out of plane hopping. In other words, even in 
the material in which the tendency to quench the orbital 
degree of freedom is strongest, there is no conduction 
band mediated antiferromagnetism in the physically rel- 
evant regime. We believe, therefore, that the physics of 
n « Uorb is simply not relevant to the manganites, in 
contrast to the assertions made iriCI. 

In j±he same way, the phase separation discussed in 
Refs.0 is crucially dependent on the existence of a com- 
mensurate antiferromagnetic order; the (0, 0, tt) order ob- 
served in the actual materials would not lead to the same 
sort of phase separation, because it would permit metallic 
in-plane conduction. 
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B. Numerical results 

By solving the dynamical mean field equations numer- 
ically, we have computed the ferromagnetic transition 
temperature Tc as a function of the Hund's coupling. 
Fig. 1^ shows ferromagnetic vs J for fillings of n = 0.7 
(solid circles) and n = 0.25 (squares), for doubly degen- 
erate eg orbitals {uorb = '2;n/norb = 0.35 or 0.125). As 
we noted in the analytic treatment of the low J limit, 
for not too low densities {n > 0.2norb), a ferromagnetic 
solution cannot be sustained. Thus for a modest n, as 
J is decreased the ferromagnetic Tc vanishes and the 
ground state changes from ferromagnetic to incommen- 
surate to antiferromagnetic. As noted previously, we find 
Uc ~ 0.195no,b ~ 0.2norb- We therefore suspect that tbe 
n = 0.2{norb = 1) curve in the Tc — J diagram of Ref.H, 
which indicates that Tc drops to at a finite J of order 
1, is incorrect at the low J end; the FM solution should 
be sustainable down to a very small J. 
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FIG. 4. Tc vs J plot of the double-exchange model for 
doubly degenerate Eg orbitals. The n = 0.7; Uorb = 2 curve 
is relevant for Lao.yAo.sMnOs compounds. For n < 0.39, the 
model can sustain a ferromagnetic solution down to the lowest 
J. The inset amplifies the low J end of the n = 0.25 curve, 
to illustrate how the numerics (open boxes) patch on to the 
analytical expression (stars) from Eq. (^4|). The solid lines 
are a guide to the eye. Tc{J — oo) = 0.079 for n — 0.7 and 
0.055 for n = 0.25. 

The stars in Fig. ^ correspond to Eq. (p4|), with the 
chemical potential fj, = — 1.52i for n = 0.25;n/norb = 
0.125(the curve is magnified in the inset). We notice that 
one cannot sustain a ferromagnetic transition for n = 0.7 
as J ^ 0. 

The ferromagnetic transition temperature has an ap- 
parently complicated dependence on interaction strength 
and doping. Wc now show that this seemingly compli- 
cated dependence simply reflects the variation of the ki- 



netic energy, K, with these parameters. Fig. || plots Tc 
against the change in kinetic energy between the para- 
magnetic and the T = ferromagnetic state {AK = 
K{Q) — K{Tc)), both being normalized by the K of the 
noninteracting (J = 0) system at T = 0. 
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FIG. 5. AK = K{0) - K{Tc) vs. n for the DE model for 
n — 0.7; Uorb = 2. Both the quantities are scaled with the K 
for the non-interacting system at T = (A'j=o(T = 0) = l.Olt 
for n = 0.7). 

To understand this curve, consider the general expres- 
sion for the free energy, F ~ E - TS of -ff_D_B (Eq. |^), 
where S is the entropy and E — —K — J where K is 
given by Eq. (H) and 

J ^S^■d\^ac.^id^^) (28) 

= [6„^cr - Jo-] . 

n 

where the latter expression is the dynamical mean- 
field result. Entropy favours the disordered state; fer- 
romagnetic order is driven by a decrease in energy. 
Tc > implies i?/m < -E-nm, where -E/m is the 
ground state energy of the ferromagnet and Enm that of 
" non- ferromagnet" (incommensurate or antiferromagetic 
order). At low T all spins are aligned, and carriers can 
move freely, while at high T the random spin arrange- 
ment and large J means carrier hopping is somewhat 
blocked. Tc is therefore set by the change in K(T^, which 
is a simple number (« 1/3) of K{T = 0) as J oo. 

As J ^ oo, the local spin is always parallel to the core 
spin and J does not change between ferromagnetic and 
paramagnetic states. The transition is entirely driven by 
the change in K between the paramagnetic and ferro- 
magnetic states. In the non-ferromagnetic (NM) phase 
the band is narrower than in the FM phase, but at fi- 
nite J the electron has some possibility to hop onto the 
"wrong" spin site. This is equivalent to the on-site mag- 



7 



netization not being saturated in the NM phase, thus 

j7/m > Jnm- SluCe Tc ~ E^m - Efm = AK + AJ, 

it impHes that Efm = Enm when Tc = and we get 
Kfrn < Knm at this point. So at high J, we start with 
fully polarized bands and with the magnetic state having 
more kinetic energy. The high J expansion suggests that 
everything comes from virtual hops : AK ^ / J and 
AJ ~ J An ^ J ■ t/ J ^ t. Since AJ is independent of J 
at large J,Tc is linear in AK for large values of J. As we 
reduce J to a value where the polarization isn't complete, 
first K^s cross and then Tj, — > 0. We have used the dy- 
namical mean-field equations to determine the points at 
which Efm = Enm iu Fig. ||. For n — 0.7{n/norb — 0.35) 
the critical J O.Qt. 

Since K is of crucial interest to us, we present its val- 
ues at T = for the soluble limits J = 0, oo in Table I 
for a variety of fillings. At J = the system is always 
paramagnetic, and the kinetic energy increases with n in 
the range < n < 1. For J — oo however, the kinetic 
energy is maximal for n = Uorb/'^- As K(Tc) — K{0)/^/2 
for J = oo and the transition temperature is tied to 
the kinetic energy (Fig. also is a maximurB-jajt 

the same filling, as has been noted in earlier workstZlEl. 



TABLE I. The kinetic energy K a,t T = 0, evaluated using 
a dynamical mean field method with a semicircular density of 
states, for J — and J = oo at different fillings. 



n/norb 


A'( J = 0)/n„rt 


A'(J = co)/norb 


0.125 


0.2160 


0.1958 


0.25 


0.3916 


0.3248 


0.35 


0.5086 


0.3889 


0.50 


0.6496 


0.4244 


0.75 


0.7994 


0.3248 



We have argued that the kinetic energy is the crucial 
parameter determining Tc- To further substantiate this 
we show in Table II results obtained for T^LJ ^ oo) by 
a variety of techniques in a range of niodelsu, expressed 
in terms of Kj=od{T = 0). Roder et alJl3 used a series 
expansion technique. Yunoki et aLQ have studied ther- 
modynamic properties of the classical core-spin model 
in d = 3 with a single orbital nearest neighbour hop- 
ping using Monte Carlo on 6x6x6 clusters. Calderon 
and Breycl have performed similar Monte Carlo calcula- 
tions on 4x4x4 and 20x20x20 lattices, and argued that 
Yunoki et al.Q underestimated by a factor of ^ 1.6. 
The results of Calderon et al. are in agreement with 
those of Roder et al. Note that in all calculations the 
doping dependence of Tc is essentially the doping depen- 
dence of the kinetic energy. 

We now consider the dynamical mean field results. 
These yield somewhat higher Tc/K results, as expected 
because fluctuations are neglected. Calderon and Brey 
partitioned the hopping term tij into average (t) and 
random (Stij) components and argued that fluctuations 
lead to a « 25% correction to T^. Applying this cor- 



rection leads to the numbers given in brackets, which 
are in good agreement with those found in Refs.clEa. 

TABLE II. Tc/K{T = 0) for J = oo; comparison of the 
different methods. The bracketed terms for d = oo indicate Tc 
values that are ~ 25% fluctuation corrected. 



Method n/norb T^/Kj=^(T = Q) 



d=oo 0.5 0.070(0.053) 


d=™ 




0.25 0.063(0.050) 


Series 


11: 


0.5 0.053 


Series 


Li: 


0.25 0.050 


Yunoki ■ 


0.5 0.036 


Yunokit 


, 0.25 0.033 


Calderon 


I 0.5 0.056 


Calderor 


y 0.25 0.056 



VI. OPTICAL CONDUCTIVITY 

Next we discuss the optical conductivity of the double- 
exchange model. At large J the density of states of this 
model consists of two nearly semi-circular bands corre- 
sponding to conduction electrons parallel (t) and antipar- 
allel (I) to the core spin. The bands are separated by an 
energy 2 J. This structure is shown in Fig. (^, which we 

now use to give a qualitative discussion of a. 
DOS 



wrong-spin transitions 




FIG. 6. The large J density of states, consisting of bands 
of electrons parallel (|) and antiparallel (J.) to the core spin 
along with possible optical transitions indicated by arrows. 
fj, denotes the chemical potential, and the shaded area rep- 
resents the electron filling. There are two pieces to o{uj); 
the transitions within the lower band with lu ^ t and the 
"wrong-spin" transitions from the lower to the upper spin 
band around uj ~ 2J. The "wrong-spin" transitions become 
allowed as T is increased from 0. For large J, these two pieces 
are well separated. 

The crucial point is that the optical process conserves 
electron spin. At T = in the fully polarized ferromag- 
netic state all spins are aligned, making it impossible for 
an optical process to create a final state with an Cg elec- 
tron anti-aligned to a core spin. Further, the perfect spin 
alignment means that no scattering processes are present; 
thus the T = conductivity is simply given by the band 
theory of the spin polarized eg electron manifold. The 
spectral weight in these transitions follows directly from 
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the band theory kinetic energy. Impurities, weakly cou- 
pled phonons etc. will change the form of the conduc- 
tivity but will not significantly affect the total spectral 
weight. 

As T is raised, the spin disorder increases, and an elec- 
tron moved to an adjacent site may find itself anti-aligned 
to the new core spinel. 

Thus as T is raised the resulting spin disorder leads to 
a broadening of ti, a decrease in total spectral weight, and 
also to a shift of oscillator strength to the peak at ^ 2 J. 
For 2 J >> i, the wrong-spin peak is well separated from 
the parallel spin transitions and we shall derive its spec- 
tral weight in this limit. Fig. R[a) shows the real part of 
the optical conductivity (Eq.]l^), calculated within the 
dynamical mean field theory, for the experimentally rele- 
vant filling of n — 0.7, for doubly degenerate conduction 
electrons. The results shown are for J/t =1,2 and 4. 
The wrong-spin transitions produce the peak at 2 J, for 
J > 1 this becomes well separated from the same-spin 
transitions. 

We are interested here in what fraction of the total 
optical weight goes into the wrong-spin transition. The 
spectral weight ratio is insensitive to details of the band 
structure. In what follows, we give an analytical estimate 
for the oscillator strength at 2 J for large values of J and 
compare it with our numerical results. For J >> i, ^ 
2 J and T = Tc, the inter-band optical weight is 



dekN{ek)(p{ek) I — A^{ek,uj) 



TT J 



(29) 



where we have used Eq. ( |l6| ) for a. Here, is the imag- 
inary part of the Green function for parallel electrons and 
A+ for an ti-parallel ones. , , , q ,q 
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FIG. 7. (a) The optical conductivity at T = for the 
doubly -degenerate DE model for n — 0.7. The results shown 
are for J =1,2 and 4; the inter-band peak is at 2J. (b) Spec- 
tral weight in the — Eg inter-band transitions, for n — 0.7 
and T = Tc- For J » t, the inter-band weight falls off as 
■J. The straight line has a slope of n/2, corresponding to the 
large J analytic expression Eq. (jsij). 

Now, A+ is an analytic function and its integral gives 
the real part evaluated at ui 
large J limit. Thus 



J, which is 1/(2J) in the 



de. 



(4- 



UttJ 



duj 



A^{ek,uj). 

(30) 



The w-integral gives the momentum space occupancy 
n{ek)- If the scattering is strong enough, this is just n. 
So, for large J and at the magnetic transition we get 



J ^dna2j{^) = 



n 
2J 



(31) 



which falls off as 1/J and is proportional to the filling 
n. The numerical results are shown in Fig. 0(b); we 
see that the analytic value for the slope, n/2 = 0.35 
(straight line) from Eq. 31 overlays the numerical curve 
for large values of J. This analysis shows that the 
spectral weight in the wrong-spin transitions depends 
crucially on n. A previous calculation of a used an 
norb = 1 and n = 0.7. The resulting oscillator strength 
in the wrong-spin transition is therefore not applicable 
to manganites with|-,n/norb — 0.35 (and indeed the os- 
cillator strength inEl is larger than the one we calcu- 
lated for the physical n/uorb by a factor of two). To 
aid in comparison to experiments we give in Table. Ill 
the total T = spectral weight, the total T > spec- 
tral weight and, for T > Tc, the spectral weight in the 
same-spin and wrong-spin transitions, calculated using 
the dynamical mean field theory for various values of J. 

TABLE III. The total T = spectral weight K{0), the to- 
tal T > Tc spectral weight K{Tc) and , for T > Tc, the spec- 
tral weight in the wrong-spin {Kanti)s,nd same-spin (Kpar) 
transitions for various values of J, and n = 0.35norb- 



J 


A'(0) 




Isanti 




1 


0.7778 


0.8060 


0.1209 


0.6850 


2 


0.7778 


0.6994 


0.1067 


0.5926 


3 


0.7778 


0.6541 


0.0833 


0.5707 


4 


0.7778 


0.6291 


0.0661 


0.5629 


5 


0.7778 


0.6136 


0.0547 


0.5588 


6 


0.7778 


0.6030 


0.0467 


0.5562 


7 


0.7778 


0.5952 


0.0408 


0.5543 


11 


0.7778 


0.5778 


0.0275 


0.5502 
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VII. COMPARISON WITH EXPERIMENTS 
A. Overview 

In this section we relate our results to experiments. 
The measurements we analyse are the value of the fer- 
romagnetic transition temperature and the magnitude 
and temperature dependence of the spectral weight (in- 
tegrated area) in different regions of the optical conduc- 
tivity spectrum. The physics issues we are able to clar- 
ify include the extent to which the double-exchange-only 
model describes the behaviour of the manganites, the 
proper interpretation of the optical spectrum, and the 
value of the Hund's coupling. The remainder of this sec- 
tion is organized as follows. In sub-section B wc sum- 
marize our results in a manner suited to comparison to 
data, in sub-section C we outline what is known about 
the optical conductivity data, emphasizing the heretofore 
unresolved issues arising in interpretation of experiments 
and in sub-section D we present the comparison between 
our results and data. 



B. Summary of theoretical results 

We have studied the double-exchange-only model de- 
fined by Eq. ^ This model involves itinerant electrons 
hopping among sites of a lattice and coupled ferromag- 
netically to electrically inert core spins. No other interac- 
tions are explicitly included. This model captures some 
aspects of the CMR mangapase perovskites and of other 
'double-perovskite' system^j; whether otbMj physics is 
important is the subject of present debateEZfEl. 

The double-exchange-only model was shown to have 
two crucial parameters: the electron kinetic energy K 
defined in Eq. |^ and the itinerant electron-core spin cou- 
pling J. K is the expectation value of a local oper- 
ator. It depends on parameters like temperature and 
J, but is insensitive to details of band structure. For 
the J values and carrier concentrations of physical rele- 
vance, the ground state is a fully polarized feromagnet 
and K{T = 0) is independent of J and in the double- 
exchange-only model may be computed via a simple 
band structure calculation. The value appropriate to 
Lao.rSro.sMnOa \s K ^ 0.84ey. We calculate Tc as a 
function of K{T = 0) and J using ^-jpean field method 
which treats local dynamics exactlyEj, and by compari- 
son to published Monte Carlo calculationscl were able to 
estimate the corrections to Tc from non-mean-field spa- 
tially dependent fluctuations. The dynamical mean field 
results are often published in terms of a bandwidth pa- 
ramete t; the t corresponding to the band theory K is 
tK imeV. 

By comparing our results to those obtained by other 
techniques, we found that in the J oo limit Tc/K{T = 
0) ~ 0.16, essentially independent of model details. For 
finite J, we found (for carrier concentrations relevant to 



the CMR materials with ferromagnetic ground states) 
that Tc was linearly proportional to the change, AK, 
of K between Tc and T — 0. We further found that as 
J — > oo, AK/K 1/3; as J is decreased, AK decreases. 
At a critical J, AK — and below this J a ferromagnetic 
ground state cannot be sustained. 

As previously noted, K is measurable in optical exper- 
iments, and therefore a quantitative test of the double- 
exchange-only model is possible and conversely con- 
straints on optical conductivity may be extracted from 
measured TcS. The issues involved in this analysis will be 
treated in the next sub-section. To conclude the present 
sub-section we discuss the additional information con- 
cerning the strength of non-double-exchange interactions 
which may be obtained from the comparison of Tc and 
the measured K{T). 

We may represent the electronic energy of a solid. 
Eel = —K + 1 where K was defined in Eq. and / repre- 
sents the expectation value of all interactions, including 
the J term and other interactions not included in the 
double-exchange model. In the rest of the discussion we 
assume J is large enough to have a spin-polarized ground 
state. 

In the double-exchange-only model at T = 0, if is sat- 
urated and given by the band theory value (assuming 
spin-polarized electrons). As T is increased, the entropy 
driven spin-disorder leads to a reduction in K; the change 
in K between T — and T > Tc depends on J, and be- 
comes as large as 1 — l/-\/2 ~ 30% in the J oo limit. 
Other interactions do not commute with K, and there- 
fore change the electronic state in a way which reduces 
K. Thus, we argue that K cannot be greater than the 
spin-polarized band theory value, and that a. K{T ~ 0) 
appreciably less than this value indicates other interac- 
tions are important. 

As temperature is increased from T ~ Q, spin disorder 
leads to a decrease in K and therefore to an increase in 
the relative strength of the interaction terms which in 
turn causes a further decrease in K . This self-consistent 
effect says that for a given J, the relative change AK/K 
in kinetic energy between Tc and T = increases with 
increasing intcmiition strength. This physics was inves- 
tigated in Ref.liJ for J = oo in the particular case of 
electron- phonon interaction, but is expected to be more 
general. Further, because kinetic energy is what decides 
Tc, the interaction induced decrease in K must decrease 
Tc below the value predicted by double-exchange. 

To summarize, the double-exchange-only model pre- 
dicts a definite set of relationships between Tc, K{T — 0), 
and the T-dependence of K and AK . These are summa- 
rized in Fig. H and Tables. Ill and IV. 

Adding other interactions causes Tc and K{T — 0) to 
decrease and AK to increase. 



10 



C. Optical conductivity 

This subsection discusses the interpretation of optical 
conductivity data. As we have shown, the magnitude and 
temperature dependence of the spectral weight in the eg 
contribution to the optical conductivity contains crucial 
information about the electron kinetic energy and Hund's 
coupling. In order to extract this information one must 
identify the Cg contribution to the measured conductiv- 
ity, which is not straightforward because of the overlap 
between the Cg transitions of interest and processes in- 
volving other bands. In this subsection we analyze the 
issues arising in the interpretation of the measured con- 
ductivity, with emphasis on how our theoretical results 
can be used to resolve some of the difficulties in inter- 
preting data. 

The optical conductivity cr(Lj) is the linear response 
function relating spatially uniform frequency dependent 
current i{uj) to applied electric field E{oj). In simple 
terms, the conductivity d escribes how electrons move 
in response to an electric field, and thereore contains in- 
formation about interactions which may hinder this mo- 
tion. In metals it is useful to distinguish between intra- 
conduction-band processes (those which involve scatter- 
ing of electrons between conduction band states, which 
in many cases including the CMR manganites are the 
states of immediate physical interest) and other processes 
which involve scattering of electrons from conduction 
band states to other (empty) bands, from other (filled) 
bands to the conduction bands, or which do not involve 
the conduction bands at all. The other processes are 
usually called interband, but in the manganites there are 
two orbitals per unit cell and therefore some of the intra- 
conduction-band processes are, strictly speaking, inter- 
band. 

In any event, one would like to extract from the mea- 
sured conductivity the portions pertaining to transitions 
between the states of interest (in the manganite case, the 
intra-conduction-band processes) and analyse only these. 
However in many cases involving transition metal ox- 
ides it has not been clear how to separate the interesting 
conduction band contributions from other un-interesting 
processes. The CMR materials are a promising system 
in which to investigate this issue because the conduction 
band contribution to a has a strong temperature depen- 
dence and a definite relation to Tc, which can be used to 
distinguish it from other contributions. 

The main point is this: the double-exchange-only 
model predicts a definite relation between the T = eg 
oscillator strength and Tc- Interaction corrections only 
reduce Tc below the double-exchange value and K below 
the band theory value. Therefore the part of the opti- 
cal spectrum asigned to the eg electrons must contain at 
least enough kinetic energy to reproduce the observed 
but cannot contain more kinetic energy than predicted by 
spin-polarized band theory results. 

A further constraint is provided by the change AK in 



kinetic energy between T = and T^. For Hde (Eq. |), 
AK < 0.3K{T = 0) and the maximal value occurs as 
J oo. Limited information is available concerning 
models with additional interactions, but published calcu- 
lations iar the double-exchange plus phonon problem at 
J — ^ ooEZi show that the fractional change AK/K{T = 0) 
can become slightly larger than 0.5, but the magnitude 
of AK is never much larger than the double- exchange- 
only value. This is useful because the change in optical 
spectral weight can be accurately measured, and the ob- 
served changes can with confidence be attributed to the 
Cg electrons of interest. A final constraint comes from the 
position and spectral weight of the "wrong-spin" transi- 
tions. The point is that as T ^ 0, all of the spins are 
aligned in both the ground state and all states accessi- 
ble from it via the optical matrix element. However, for 
T > Tc, the core spins are completely disordered and 
therefore when an electron hops (or is pushed via the 
optical matrix clement) from one site to another it has 
a probability for landing in the "wrong-spin" configura- 
tion, i.e. with Cg electron anti-parallel to the core. The 
probability depends on the hopping matrix element (i.e, 
the T = kinetic energy) and the value of J. This physics 
was first pointed out by Okimoto et.al3, who, however, 
obtained what we argue below was an incorrect value of 
J. We have determined the optical oscillator strength in 
this transition, as a function of K{T — 0) and J. The 
oscillator strength has a strong dependence on carrier 
density not noticed in previous work. 

D. Analysis of data 

Optical conductivity has been measji|Ujed|-iji a wide 
range of manganites by several groupsBElO'liSl. Results 
are qualitatively similar, but there are substantial quan- 
titative differences between results of different groups. 
Some of the differences seem to be eypprimental artifacts 
associated with surface preparationE2lE£l; others seem to 
relate to differences in physics among various members 
of the CMR family of materials. As we shall see, our re- 
sults provide consistency checks which allow one to sepa- 
rate artifacts from intrinsic behaviour and to make some 
statements about the underlying physics. ■— ■ 

To be concrete, we discuss the data of Quijada et.aO, 
who measured a{Lu,T) for three pseudo-cubic man- 
ganite films Lao.7Sro.3Mn03 (LSMO), Lao.TCao.sMnOa 
(LCMO) and Ndo.rSro.aMnOa (NSMO). The quahtative 
features of the data are (i) at w < 3eV, an absorption 
with a pronounced frequency dependence and an inten- 
sity that shifts to lower frequency and increases markedly 
as T is decreased from Tc to a low temperature, (ii) a 
strong feature centered at w « AeV with little apparent 
T-dependence, and (iii) a weak feature at w « 3eV, visi- 
ble as a decrease in absorption intensity as T isL/iecreased 
below Tc- The interpretations offered in Ref.tHi was that 
(1) the integral of the low T conductivity between uj = 
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and Lo = 2.7eV was a good representation of the total low 
T conduction-band spectral weight, (2) that the strong 
feature at oj ^ ieV was a Mn-0 interband transition 
and (3) the weak feature at w « 3eV was the "peak at 
J" (it is actually at 2 J in our units, as in Fig. |^) due to 
"wrong-spin" transitions characteristic of the spin disor- 
derecLstate. This identification is controversial. Okimoto 
et aLtlEl identified a lower energy feature as the peak at 
J. 

We begin our discussion with the "pi^^ at J". The 
identification offered by Okimoto et al.Ellj, implies J ~ 
0.75eV which is less than the dynamical mean field hop- 
ping t — l.OSey (this is derived from equating the dy- 
namical mean field kinetic energy K to the band theory 
K = 0.84mey). Thus the J must be larger as a ratio of 
J/t w 0.7 would imply a ferromagnetic state cannot be 
sustained at n/riorb — 0.35, in contradiction to experi- 
ment. _ 

We now turn to the identification of Quijada et al.Ej, 
which corresponds to J = 1.5eV^, i.e to J/t « 1.38 or a 
fluctuation corrected ~ 0.03t « 310iir, somewhat be- 
low the 350K observed for Lao.rSro.sMnOa. This casts 
doubt on the interpretation offered by Quijada et al. Fur- 
ther, the spectral weight observed by Quijada et al. in the 
" peak at J" , Kanti = 34meV is about twice the spectral 
weight predicted by our calculation, using the Quijada J 
and the t inferred from band theory. We conclude that 
the feature observed by Quijada et al. is unlikely to be 
the "peak at J" and that the true Hund's coupling is 
probably larger so that the peak at J is outside the ex- 
perimental range. At least one needs a value of J w 2eV 
to reproduce the observed Tc using the band theory K. 

We digress briefly on the question of the origin of the 
observed feature. It seems reasonable that it is caused 
by a shift of the energy of the Mn-0 interband transi- 
tion at 4eV, due to changes in the eg band as T is varied 
through Tc (the leading edge of this transition is an ex- 
citation from a filled oxygen state to the Fermi level in 
the e — g band). One possibility is the double-exchange 
induced shift in jphemical potential, discussed in a dif- 
ferent context ind. Unfortunately, for relevant dopings, 
the changes that we find for J/t > 2 are of the order of 
1 — 2%, too small to explain the observed feature. An al- 
ternative possibility is a polaronic shift asociated^with a 
change in effective electron-phonon interactionLjcj, but 
this has not been studied in detail. 

We now present a more detailed discussion of the re- 
mainder of the spectrum, proceeding on the assumption 
that the "peak at J" is not visible in the spectrum at 
uj < AeV, i.e J > 2eV and that the total T = spec- 
tral weight is less than or equal to the band theory value 
K = 280meV. Referring to Fig. 6 of Quijada et al.Ej, 
one observes a change between Tc and T — Q AK{uj*) — 

du>a{T ^ 0,a;) — cr(T > Tc,uj) in spectral weight of 
0.1-0.13 eV in LSMO, 0.1-0.12 eV in LCMO and 0.1 in 
NSMO. Consider LSMO first: if a double-exchange only 
model described the physics, the total T = spectral 



weight would be at least 3Aif, i.e K{T = 0) > 0.3eV. 
This is slightly greater than the band theory value. Now 
the Tc of LSMO is « 350;^ « Q.lKiT = 0). Reproduc- 
ing this Tc with K{T = 0) = Q.ieV would require a J 
slightly less than 2eV. From Table. Ill we see that such a 
J / K would imply a change AX in spectral weight much 
smaller than observed. Similar, but more severe prob- 
lems arise for the other materials, LCMO and NSMO. 
The double-exchange-only constraint K{T = 0) > AiiT 
implies K{T = 0) > O.SeV for these materials as well. 
The lower Tc's (« 270if for LCMO, « 250A' for NSMO) 
would then imply J — l.b— l.^eV , again yielding a AAT 
that is too small, and a "peak at J" visible in the spec- 
trum. 

We would like to add here that, our values of K are 
a factor of 3 higher than that oitlj. This is because of 
the l/Vd dependence of each current vertex in the ex- 
pression for (j{u>). Thus, for d = 3, the kinetic energies 
of Quijada et al should be multiplied by (\/d)^ = 3 to 
compare with our numbers. We do this in Table |l^ for 
facilitate comparison. 

To conclude this section we invert the logic, using 
the observed transition temperatures and AAT to obtain 
bounds and estimates for the conduction band spectral 
weights. Wc argue that for a given K{T = 0) and J, 
the double-exchange-only model gives an upper bound 
for Tc. Thus values of Tc and J yield a lower bound for 
K{T = 0). These bounds are shown in the Table IV for 
J = 2eV , ieV and AeV for the systems studied by Qui- 
jada et al. Note that an upper bound of K{T = 0) is 
given by the band theory value. 

An alternative estimate may be obtained by con- 
sidering the limit J ^ oo. In this case, as noted 
inE3, Tc is a universal function of K(Tc) (which may 
itself be affected by other interactions). Further, Tc 
must decrease as J is decreased; thus, we may ob- 
tain a lower bound on K{Tc) from the J = oo re- 
sult for Tc- These bounds are also listed in Table IV. 



TABLE IV. The lower bounds on the kinetic energies 
(Kmin) at T = and Tc for the systems studied by Quijada 
et al., shown for J = 2eV,3eV and ieV. The upper bound 
is given by the band theory value Kband{0) = 840?i*«y. To 
compare with experiments, we cite the values of Ref.LJ in the 
last two rows. The experimental values for K arp*epresented 
in the conventions of this paper, i.e. values fromli2l multiplied 
by 3. All energies are in meV. 





LSMO 


LCMO 


NSMO 


Krmn{0); J = 2eV 


630 


392 


350 


A"m«„(0); J = 3eV 


535 


336 


327 


A'™„(0); J = 4eV 


514 


345 


314 


Kmin(Tc); J = oo 


477 


330 


309 


Kexpt (0) 


780 


660 


600 


t^expt (Tc) 


477 


393 


390 



It is interesting to note that the experimental K{Tc) 
for LSMO is 477meV, which exactly saturates the lower 
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bound on K{Tc) from double- exchange. Because the 
infinite- J double-exchange-only is expected to be the 
upper bound to the true Tc of an interacting model we 
expect that in the actual material the J is larger than 
2eV (so the model is not far from the J = oo limit and 
we suspect that the integration up to = 2.7eV does 
not capture quite all of the spectral weight. 



VIII. CONCLUSIONS 

We have given a complete and correct treatment of the 
phase diagram, kinetic energy and some aspects of the 
optical conductivity of the doublc-exchange-only model 
of mobile carriers coupled to core spins. We have deter- 
mined the physics operating in different regions of the 
phase diagram and have demonstrated the importance 
of choosing parameters (especially carrier density) ap- 
propriate to the materials of interest by exhibiting the 
incorrect results obtained by the use of wrongly chosen 
parameters. We have shown that the crucial quantity is 
the electron kinetic energy and have used our results to 
estimate the kinetic energy of several manganite systems. 
Our results also provide insight into the crucial question 
of which portions of the spectrum pertain to the low en- 
ergy electronic degrees of freedom. 

There are several directions for future work. One is 
to combine the dynamical mean field method with a re- 
alistic band structure, to obtain a treatment of the fre- 
quency dependence of a. Another is to employ the meth- 
ods presented here to models where double-exchange is 
combined with other interactions. If this were carried 
through, it seems likely that the combination of Tc and 
the changes in optical spectral weight could be analysed 
to provide detailed knowledge of the strength, energy 
scale, and nature of any additional couplings. 
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